Topology of amorphous tetrahedral semiconductors on intermediate lengthscales 
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Using the recently-proposed "activation-relaxation technique" for optimizing complex structures, 
we develop a structural model appropriate to a-GaAs which is almost free of odd-membered rings, 
i.e., wrong bonds, and possesses an almost perfect coordination of four. The model is found to 
be superior to structures obtained from much more computer-intensive tight-binding or quantum 
molecular-dynamics simulations. For the elemental system a-Si, where wrong bonds do not exist, 
the cost in elastic energy for removing odd-membered rings is such that the traditional continuous- 
random network is appropriate. Our study thus provides, for the first time, direct information on 
the nature of intermediate-range topology in amorphous tetrahedral semiconductors. 



The structure of amorphous materials on short and 
intermediate lengthscales remains largely, unresolved in 
spite of decades of work on the problemJjlJ While diffrac- 
tion experiments can in principle provide the desired in- 
formation, they still lack the sensitivity necessary to dis- 
tinguish between various possible conformations of the 
structure. They must be backed up, therefore, by accu- 
rate structural models in order to allow a meaningful in- 
terpretation of the data. Models, however, are difficult to 
construct because they require an accurate description of 
the interatomic forces, which in turn limits the timescale 
over which simulations can be carried out. Since glasses 
relax very slowly, this is an extremely serious problem 
that can only be addressed through the use of judicious 
optimization methods. 

Amorphous tetrahedral semiconductors, either elemen- 
tal or compound, are potential candidates for the fabri- 
cation of micro- and opto-electronic devices. Advances 
in this area have however been hindered by the lim- 
ited understanding of their structure, and thus electronic 
properties, even at the most elementary level. For in- 
stance, in the case of amorphous silicon (a-Si), the aver- 
age number of nearest-neighbors is still not known pre- 
cisely, even though evidence (experimental and theoret- 
ical) accumulates that it is very close to 4. Thus, it 
would appear that the fundamental "building block" of 
these materials is a tetrahedron — a single atom and its 
four nearest-neighbors arranged in a pattern that closely 
resembles that found in the crystals. It is in the ar- 
rangement of those tetrahedra, i.e., on intermediate (and 
larger) lengthscales, that differences between amorphous 
and crystalline materials arisa§B and on which, therefore, 
attention should be focused. 

The case of compound III-V materials, such as GaAs 
and InP, is more problematic because the building blocks 
are not as clearly defined. Indeed, the chemical iden- 
tity of the atoms constituting the tetrahedron may vary: 
While in the zinc-blende crystal each atom of a given 
species is surrounded by four atoms of the other species, 



this is not necessarily so in the amorphous material, be- 
cause of disorder. Thus, there is a possibility that "wrong 
bonds" be present in the amorphous phase, i.e., bonds 
between like atoms. This will result in some distortion of 
the tetrahedra, and therefore on their relative arrange- 
ment. However, III-V compounds are partly ionic and 
the cost in energy of wrong bonds is expected to be sig- 
nificant; the number of them, therefore, will be deter- 
mined by the competition between elastic deformation 
and chemical disorder. There is currently no expeximcn- 
tal method for measuring this number in a-GaAs.u 

In this Letter, we aim at determining the optimal den- 
sity of wrong bonds in amorphous III-V semiconductors 
through a comparison with elemental materials, where 
such defects do not exist. We illustrate our ideas by con- 
sidering the prototypical materials a-Si and a-GaAs. We 
note first that wrong bonds are closely related to the 
presence of "odd-membered rings" in the structure, i.e., 
closed paths between an atom and itself containing an 
odd number of bonds. Evidently, only even-membered 
rings are present in zinc-blende crystals. Also, it is al- 
ways possible to "decorate" an amorphous network con- 
taining only even rings with two types of atoms such 
that there are no wrong bonds. Thus, we may formulate 
the question in the following manner: are odd-membered 
rings energetically favorable, i.e., necessary to the sta- 
bility of amorphous tetrahedral semiconductors, be they 
elemental or compound? And is it possible to construct 
a tetrahedrally-bonded amorphous model without odd- 
membered rings? 

Anticipating our results, we find that it is indeed quite 
possible, albeit difficult, to construct such an all-even- 
ring network, but the cost in elastic-deformation energy 
required by this constraint is such that the resulting 
model is not appropriate to elemental materials. In con- 
trast, a model that contains odd rings is found to be 
unappropriate for III-V materials, because the cost in en- 
ergy associated with wrong bonds largely exceeds that of 
the elastic deformation. Thus, odd-membered rings must 
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be present in a-Si and absent in a-GaAs. This consti- 
tutes, to our knowledge, the first evidence of a topological 
difference (i.e., without regards to chemical identity) be- 
tween the two materials. We find, also, that the presence 
(or absence) of odd-membered rings is extremely difficult 
to extract from such quantities as the radial distribution 
function or distribution of dihedral angles, and therefore 
realistic structural models are a necessity for proper in- 
terpretation of experimental data. 

In order to establish these results, we proceed as fol- 
lows. Using the "activation-relaxation" technique^ (see 
below), we construct two "generic" continuous-random 
networks (CRN)El with periodic boundary conditions: 
The first, CRN- A, .possesses the same attributes as an 
infinite Polk modelo as constructed by Wooten, Winer 
and Weairep and thus contains odd-membered rings. (In 
fact, the WWW algorithm consists in amorphizing a crys- 
talline lattice through, precisely, the introduction of odd 
rings). Our second model, CRN-B, was constructed with 
the constraint that odd rings be absent, and thus corre- 
sponds to an infinite Connell-Temkiro computer model. 
(The Connell-Temkin model is a hand-built, 238-atom, 
mechanical model with free surfaces). Both CRN- A and 
CRN-B were then "decorated" with either Si or Ga and 
As — so that we have, altogether, 4 different struc- 
tural models — and relaxed to their equilibrium state 
using conjugate gradients with forces derived from tight- 
binding (TB) potentials. 

It is important to note that there have been numer- 
ous attempts at simulating the structure of a-Si using 
molecular dynamics (MD) and a variety of interatomic 
potentials — from empirical to fully ab initio. All lead to 
a sizable proportion of odd-membered rings. For GaAs, 
for which there exists no satisfactory empirical potential, 
calculations have been carried out using either TBE30 or 
Car-Parrinelldl3 (CP) MD. Again, here, odd-membered 
rings — and therefore wrong bonds — are present, in an 
amount which appears to be quite large (10-12%). Since 
the period of time that can be simulated using such ap- 
proaches is inevitably short (of the order of 10-100 ps), 
there remains a possibility that the large density of de- 
fects is the result of incomplete relaxation. In order to 
resolve this point, a different, much more efficient, opti- 
mization method must therefore be used to prepare the 
structural models. 

The activation-relaxation Monte- 

Carlo technique. (ART), recently proposed by Barkema 
and Mousseaujj, provides a powerful and efficient way 
of searching configuration space for a global minimum. 
Full details of the method can be found in Ref. ^; it 
can be summarized as follows: Starting with a random 
distribution of atoms, the system is first relaxed to the 
nearest local minimum on the potential energy surface 
using standard minimization techniques. It is then "acti- 
vated" by pushing it to a neighboring saddle point along 
a valley, i.e., a path of minimal energy. The system is 
then relaxed, again, to the nearest local minimum, and 
the process iterated until a global minimum is found. 



The potential energy for the ART preparatiotupf CRN- 
A was taken to be of the Stillinger- Weber formt^l, appro- 
priate to Si, except that the strength of the bond-bending 
term was increased by 50% so as to provide a better de- 
scription of the materialE For CRN-B, the same poten- 
tial was used — the lattice parameters of c-Si and c-GaAs 
are nearly the same — but a repulsive term between like 
atoms was added in order to describe the chemical order: 
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where the sum is over all pairs of atoms, e is the Stillingcr- 
Weber energy parameter, Aij — 1.2 for like particles and 
zero otherwise, and = 3.6 A. Both CRN-A and CRN- 
B models contain 216 atoms in a cubic cell with periodic 
boundary conditions. The ART relaxation took about 
three days on a RS8000 workstation for CRN-A and 
about a week for CRN-B. This should be contrasted with 
several weeks for a TB simulation of a 64-atom system 
on a workstation or a CP simulation on a state-of-the-art 
parallel supercomputer. 

After reaching a stable configuration in the ART sim- 
ulations, the two models were made material-specific by 
further relaxing under TB potentials. Both CRN-A and 
CRN-B models wereuthus relaxed with both the Vogl et al. 
TB potential for Siti-and the Molteni-Colombo-Miglio 
TB model for GaAsB In the case of model CRN-A- 
GaAs, it is necessary to "label" the atoms. This was done 
by assigning identities at random, then iteratively mini- 
mizing the number of wrong bonds by a "label-switching" 
process. This lead to a proportion of wrong bonds of ap- 
proximately 14%, corresponding to the theoretical value 
for optimal ordering on a Polk-type network.li3. 

The energies of our four different TB-relaxed models 
are listed in Table Q. It was observed that the topology 
of the models is little affected by the final TB relaxation, 
indicating that the ART simulations effectively took the 
systems very close to the global minimum. 

Considering Si first, it is clear that the CRN-A model, 
which does contain odd-membered rings, is a much more 
favorable structure for the material than CRN-B, which 
has no such rings. The cost in elastic energy required 
to eliminate odd-membered rings can be extracted from 
Table Q as the difference in energy between CRN-A and 
CRN-B for Si, and comes out to be about 0.08 eV/atom. 

For GaAs, now, we find the opposite result: As indi- 
cated in Table |, it is much more interesting for this 

TABLE I. Energy per atom of the two networks relaxed 
with TB. For GaAs, we also give the results from the TB-MD 
simulations of Seong and Lewis (SL), Ref. [ll]. 



Network 


TB parameters 






Si 


GaAs 


GaAs (SL) 


CRN-A 


-12.735 


-13.450 




CRN-B 


-12.659 


-13.561 


-13.450 


Crystal 


-13.084 


-13.802 


-13.802 
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material to adopt the CRN-B structure, free of odd- 
membered rings and therefore wrong bonds, than the 
CRN-A structure. Thus, the energy gained by the sup- 
pression of wrong bonds largely exceeds the cost for elasti- 
cally deforming the network to eliminate odd-membered 
rings — 0.08 ev/atom, as we have just seen. 

The structural properties of the two optimal models — 
CRN-A for Si and CRN-B for GaAs — are summarized 
in Table ||. Evidently, for Si, we achieve an excellent- 
quality model, with an almost perfect coordination of 4, 
sharply peaked at this value, i.e., a structure which, while 
amorphous (see below), is almost free of topological de- 
fects. 

Likewise, our CRN-B-GaAs model, with less than 4% 
of wrong bonds and an almost perfect coordination of 
4, is by far the closest realization of a perfect, infinite, 
Connell-Temkin-like CRN ever achieved. The present 
model compares very iavorably with experimental static- 
structure-factor datajfl as demonstrated in Fig. [[] as well 
as with the results from TB- or CP-MD simulations, as 
can be seen in Tables | and || Using MD and exactly 
the same TB model, Seong and LewisO (SL) obtained an 
amorphous model whose energy lies significantly above 
that of the present modeL-, (It si not possible to com- 
pare energies from the CPt3 and TB models). The den- 
sity of wrong bonds in the present model is also much 
lower than that obtained in the TB- and CP-MD sim- 
ulations. In addition, the distribution of coordination 
numbers here is sharply peaked at about 4, while it is 
much broader in the other models. These results sug- 
gest that the usual melt- and- quench MD approach may 
not yield a fully-relaxed configuration, while ART simu- 
lation approaches this point much more effectively. ART 
relaxation under TB or first-principles forces is however 
presently not available. 

It is of interest to examine our models in terms of 
a possible structural signature of the presence of odd- 
membered rings (or wrong-bonds), not based on knowl- 
edge of the system at the atomic level, that would pro- 
vide (e.g., from experiment) information on intermediate- 
range correlations. One possibility is the radial distri- 

TABLE II. Distribution of coordination numbers Z and 
average value, nearest-neighbor distance rjvjv, and density of 
wrong bonds (WB) in our two optimized models, compared 
with other models — SL: TB-MD of Ref. [j]} CP: CP-MD of 
Ref. |l§ For CRN-A, the density of wrong bonds is that of 
the GaAs-decorated lattice as discussed in the text. 



CRN-A-Si CRN-B-GaAs GaAs (SL) GaAs (CP) 



Z = 3 


0.037 


0.051 


0.242 


0.219 


Z = 4 


0.963 


0.944 


0.598 


0.781 


Z = 5 





0.005 


0.129 





Z = 6 








0.024 





Z = 7 








0.007 
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FIG. 1. Total static structure factors for the CRN-B-GaAs 
model (full line) and from experiment (dots, Ref. ^|). 

bution function, G(r), which we show in Fig. || for our 
two best models, CRN-A-Si and CRN-B-GaAs. In spite 
of the fact that these models differ markedly in their 
medium-range topology, as we have seen, the differences 
between the two curves are clearly very small. This is a 
bit of a surprise in that it contradicts the current view 
that differences between the two structures should be 
clearly visible beyond the second peakJl3 In fact, it was 
argued by Connell and TemkinB that the signature of a 
CRN-B-type network would be, in comparison to CRN- 
A, (i) a broader second peak, (ii) a deeper minimum be- 
tween the 2nd and 3rd peak as well as (iii) a shifted and 
(iv) broader 3rd peak. This is not what we find here, 
as can be seen in Fig. ||: although there are some dif- 
ferences — not necessarily in agreement with the above 
criteria — they remain small and could easily be within 
the accuracy of the models. Thus, G(r) exhibits no clear 
signature of the presence of odd-membered rings (even 
though the differences we observe between the two struc- 
tures are within experimental resolution). 

We show, in the inset to Fig. ||, the distribution of 
dihedral angles (i.e., between second nearest-neighbc* 
bonds) in our two models. In the Connell- Temkin modelu 
(corresponding to our CRN-B), staggered configurations 
(4> = 60°) are found ta be much more probable than 
in the usual Polk CRN,cl which contains odd-membered 
rings. Our models demonstrate that the differences are 
in fact very small and again not significant. In particular 
both models seem to like staggered configurations with 
comparable probability. Thus, the distribution of dihe- 
dral angles is not a good measure of intermediate-range 
correlations. 

It follows from the above discussion that the nature 
of intermediate-range correlations cannot easily be de- 
termined with current experimental techniques, so that 
atomic-level models, which have the power to resolve dif- 
ferent possible structures yielding the same structure 
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FIG. 2. Radial distribution functions for the CRN-A-Si 
(dashed line) and the CRN-B-GaAs (full line) models. The 
GaAs curve has been rescaled along r so as to match the 
first-neighbor peak of Si. The inset shows the distribution of 
dihedral angles for the same two models. 

factor on the basis of total-energy calculations, remain 
the only resort for proper interpretation of experimen- 
tal data. Thus, while such a method as "reverse Monte 
Carlo" (Ref. |l7|) can provide almost perfect agreement 
with diffraction data, it clearly cannot resolve two struc- 
tures that yield comparable structure factors because it 
does not include a total-energy component. 

In conclusion, using the newly-proposed activation- 
relaxation technique, we have been able to construct a 
structural model appropriate to a-GaAs that contains 
almost no wrong bonds (i.e., is free of odd-membered 
rings) with an almost perfect coordination of four. This 
Connell-Temkin-like model is found to be superior, from 
both structural and energetic viewpoints, to structures 
obtained from TB and quantum MD simulations, at a 
fraction of the computational cost. In contrast, in the 
elemental system a-Si, where wrong bonds do not exist, 
the absence of odd-membered rings costs elastic energy, 
so that the traditional Polk-like CRN is appropriate. Our 
study thus provides, for the first time, direct and solid 
information on the nature of intermediate-range correla- 
tions in amorphous tetrahedral semiconductors, and, in 
particular, establishes clearly that these materials cannot 
be described by a single topological model. It underlines, 
moreover, the fact that state-of-the-art total-energy com- 
puter models are essential for reliable interpretation of 
experimental diffraction (and other) data. Finally, our 
study gives unambiguous evidence for the ability of ART 
to yield the optimal structure of strongly-disordered sys- 
tems, thus opening a new way into "complexity" . 
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